Creation of macroscopic superposition states from arrays of Bose-Einstein condensates 
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We consider how macroscopic quantum superpositions may be created from arrays of Bose- 
Einstein condensates. We study a system of three condensates in Fock states, all with the same 
number of atoms and show that this has the form of a highly entangled superposition of different 
quasi-momenta. We then show how, by partially releasing these condensates and detecting an inter- 
ference pattern where they overlap, it is possible to create a macroscopic superposition of different 

, relative phases for the remaining portions of the condensates. We discuss methods for confirming 

■ these superpositions. 

o : 

^ , PACS numbers: 03.75.Gg, 03.65.Ta, 03.75.Lm 

^ ' Quantum mechanics allows objects to exist in a coherent superposition of different states. This does not depend 
on the size of the system and means that it should be possible to create superpositions of macroscopically distinct 
states. Schrodinger first put forward this idea in his well-known thought experiment entangling the fate of a cat with 
the state of a radioactively decaying nucleus H . In this system, the cat is in a superposition of being alive and dead 
depending on the state of the nucleus. This outcome is at odds with our classical perception of the physical world 
T-H ' and has generated a great deal of discussion about whether macroscopic superposition states can be observed. 
K*" , Under conditions that carefully avoid decoherence , it is possible to create superposition states in the laboratory. 
' As time has gone on, these superpositions have been demonstrated in systems of increasing size including the spatial 
, coordinates of a single atom 0] , the quantum phase of optical coherent states |^ , the position of fuUerene molecules 
' Q , and the current in a superconductor . 

I As we shall see, Bose-Einstein condensates are a promising starting point for creating and studying cat states since 
\^ ' they consist of a large number of particles all in the same quantum state and are sufficiently cold to enable quantum 
' phase transitions. They can be considered to be analogous to lasers, which are useful for creating non-classical states 
, of light. There have been a number of theoretical proposals for how macroscopic superpositions of condensates may 
be achieved 7]. In this paper, we focus on a system of condensates transferred to an optical lattice H,I3- Although 
^ ' , , some authors reserve the term 'cat state' for a macroscopic superposition entangled with a microscopic state in the 
(~| ■ sense of Schrodinger's original thought experiment, we will use the term more generally in this paper to mean any 
Ci ' superposition of macroscopically distinct states. 

As our starting point we consider states where each site in the lattice has precisely the same number of atoms. 
We begin by showing that this state can be considered to be cat-like when viewed in quasi-momentum space. We 
then show how, by measuring interference patterns between overlapping clouds of this state, it is possible to create 
macroscopic superpositions of different relative phases between the condensates. The initial state may be created by 
^ I , a Mott insulator transition and it is helpful to begin by reviewing this transition for atomic condensates. 

An array of condensates [lo| can be created by applying, to a single trapped condensate, an optical lattice in the 
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form of a standing wave. If this is done adiabatically in the sense that we remain in the ground state, the system can 
be described by the Bose- Hubbard Hamiltonian 

<i,j> i 

where Oj is the annihilation operator for an atom at site i and, in the first term, the summation is taken over nearest 
neighbors. The strength of the tunneling between sites, J, can be adjusted in experiments by changing the intensity 
of the standing wave, thus altering the potential barrier between sites. The interaction strength between atoms, U, 
is at best only weakly dependent on the potential, but can be controlled by usin g Fe shbach resonances. As U/J 
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is increased, the atom number fluctuations at each site are progressively reduced as has been experimentally 
demonstrated M- For sufflciently large values of C// J, a quantum phase transition to the Mott insulator state has 
been predicted and observed In the limit U/J — > oo, each site has precisely the same number of atoms 
assuming commensurability, i.e. the ratio of atoms to sites is an integer. 

We shall take this perfectly-squeezed n-tuple Fock state, where n is the number of sites, as our starting point. 
The simplest system required for our scheme is n = 3 and we shall consider this case in detail. Three-well arrays of 
condensates have already been demonstrated in the laboratory Later in the paper, we will consider n = 4 and 
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comment on the applicability of our results to larger arrays of condensates. For n = 3, we label the atom annihilation 
operators for each condensate a, b, and c, respectively and take each condensate to initially have N atoms. The state 
of the system can then be written in the atom number basis as, 

\^) = \N,N,N)abc. (2) 

This is an idealized version of the state created in Mott transition experiments and we will consider below the effect 
of deviations from this due to imperfect squeezing. 

To begin with, we would like to consider the form of this state ^ when it is transformed to a quasi-momentum 
basis, i.e. a complete basis of states that have equal weightings of a, b, and c and a linearly varying phase across 
them. The annihilation operators, {aj,/?j,7^}, corresponding to a general orthonormal quasi-momentum basis for 
three lattice sites are 

-ib + e-*2«c) (3) 

-»(4-27r/3)^_^g-j(25-47r/3)^^ (-^^ 

where the angle ^ can take any value. This basis also corresponds to the outputs from six-port beam splitters ('fritters') 
if a, b, and c are the inputs These have been experimentally realized for photons [Ij]. 
In this new basis, the initial state 1-0) is given by, 

IV') - -^{a^b^C^f\0)atc 



4' + + 7f - Sal^l^l |O)„0^, (6) 
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where \0)abc and |0)c(^^ denote the vacuum states in the two bases. The state © can be written in the form 

3JV 3N-m 



1-0)=^ ^ f{m,n)\m,n,3N ~m~n)ai3j, (7) 

m— n— 

where /(m, n) are the coefficients found by operating on the vacuum state in ©. The probability of finding Na atoms 
in mode a, N/^ atoms in mode /3 and hence 3A^ — Na — N/s atoms in mode 7, is then 

P{Na,Np)^\fiNa,N^)\\ (8) 

This probability distribution is plotted in Figure 1 for 3N = 60 and is independent of the value of ^ chosen in 
the mode decomposition. This state is analogous to that formed by passing Fock state pairs through a 50:50 beam 
splitter It could be considered to have a 'Schrodinger cat-like' form in that the atoms tend to cluster around the 
corners, which correspond to all the atoms having the same quasi-momentum. A 'true cat' in the sense that we use in 
this paper would consist of a coherent superposition of all the atoms being in a and all in /3 and all in 7. In the other 
extreme, if our state was \Tp) — \N, 0, 0)abc, this would give (x (a^ + /3| -f 7^)^|0, 0, 0)af3-f in the quasi-momentum 
basis, i.e. it would be a microscopic superposition of different quasi-momenta. The state given by lO and depicted 
in Figure 1 falls between these two extremes and we will call it a three-cornered 'hat state'. 

So far we have studied only the pure number states that emerge in the limit of strong repulsive interactions and 
weak tunneling. It is interesting to consider the effect of varying the interaction strength. By analyzing the ground 
state that results from numerically diagonalizing the Bose-Hubbard Hamiltonian J^l, we can investigate how the 
distinct hat state structure evolves as a function of U/J. 

Figure 2(a) shows the probability distribution © for the ground state of the Bose-Hubbard Hamiltonian with 
U/J = 50. The interactions lead to a depiction of the condensate mode aj, i.e. the ground state contains admixtures 
of states with a few particles in modes /?{ and 75. If the ground state were a single pure condensate, as for U/J = 0, 
Figure 2 would show a single peak at (NcN/s) — (60,0). As one approaches the n-tuple Fock state regime, the 
distinct structure of Figure 1 starts to emerge. For U/J = 500, for example, the superposition of macroscopically 
occupied modes becomes clear (see Figure 2(b)). The symmetric three-cornered hat structure emerges only for the 
triple Fock state, which is to say in the limit of highly number squeezed systems. 
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FIG. 1: The probability distribution given by Q for the number of atoms in the quasi-momentum modes a and P for a state 
of the form |i/>) = |iV, N, N)abc where 3iV = 60. 



States very close to our n-tuple Fock states have been experimentally created in large lattices |9j and it would be 
interesting to see how they could be used. They may be well-suited to high precision measurement schemes in analogy 
with their two-site counterparts |1?tI In order to optimize proposed measurement schemes, we would like to start 
with a system with a large number of atoms in each lattice site since the phase resolution scales as l/{dy^), where n 
is the number of sites and d is the number of particles per site ^3 • It is difficult to achieve the Mott transition in this 
regime |l8j| . One possibility may be to first perform the Mott transition in three dimensions with only a small number 
of atoms per site. The confining potential in two orthogonal directions could then be lowered sufficiently slowly that 
the system remains in the ground state of the evolving Hamiltonian. This would leave us with a one-dimensional array 
of condensates with the same (large) number of atoms in each site if the lattice had a cubic symmetry. The influence 
of different symmetries is an interesting question that would need to be addressed in any experimental implementation 
of this scheme. 

We have seen how a cat-like state emerges in the quasi-momentum basis as the ratio U/J is increased. Now we 
would like to turn our attention to how a 'true' cat state can be created from this resource. For simplicity, we will 
consider only the case that we begin with a perfectly squeezed state. We shall show that, by partially releasing a 
triple Fock state Q and detecting an interference pattern between the expanding atomic clouds, a cat state can be 
created in the relative phase of the remaining condensates. This is an extension of work that showed that a single 
relative phase is created between a pair of condensates when atoms are detected where they overlap 0, |23, 0| . The 




FIG. 2; Probability distribution JHJ for the ground state of the Bose-Hubbard Hamiltonian with 3 sites and 60 particles for (a) 
U/J = 50 and (b) U/J = 500. 
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FIG. 3: General scheme for generating cats in relative phase space by observing interference between three condensates. In the 
far field Ai? > 



system we consider consists of three trapped condensates each initially containing N atoms and spaced from each 
other by a distance, d. At time t — 0, the interactions are switched off and the condensates are partially released from 
the trapping potential and allowed to expand and overlap. This release could be achieved, for example, by coherently 
transferring part of the population to an untrapped state. The overlapping atoms interfere and are detected at a 
distance R from the lattice (see Figure 3). 

The annihilation operator, 17, that corresponds to detecting an atom at angle 9 and distance i? ^ d at time t > 

is, 

n^^(a e-^^^'l^ + h e-'-S'-Vfi + ^ e'^^^*/'') , (9) 

where Ea, and Ec are respectively the energies required for atoms from a, b, and c to have the correct speed to 
reach the screen at 9 in time t. These kinetic energies, Ej, correspond to de Broglie wavelengths, Aj, and can be 
rewritten as, Ej — m_R|/(2i^) for j G {a, 6, c}, where Ra, Rb, and Rc are the distances travelled by an atom from 
a, b, and c to the point of detection and m is the mass of an atom. For angles such that sinfl 3> d/R, we can write 
{Ra, Rb, Rc) — (R ~ ds'm9, R, R + dsm9). However, we need to ensure that we detect atoms at sufficiently small 
angles to correspond to the scale of variation for the interference fringes, i.e. sm9 ^ A/d, where A is the typical de 
Broglie wavelength of the atoms. This gives us a condition for the geometry of the set-up, X/d^ d/R, i.e. XR ^ d^, 
which is the usual far-field condition guaranteeing that the size of the diffraction pattern is large compared with the 
source size. Substituting into @, we obtain, 

17 ^ _L (a + 5 e""' + ce-^*""'"") , (10) 

where rj — mRd/Tit and we have ignored any overall phase. Since rj depends on t, if many particles are detected as 
they arrive at a screen at various times, we get a range of values of 77 and the fringes will tend to wash out. To avoid 
this, we can take a 'snap-shot' of the spatial distribution of the overlapping clouds at a given time. The time we need 
to wait in order to satisfy the far- field condition is t 3> md^ /h. 

For simplicity, we will consider that atoms are detected at small angles and so the replacement sin 9^9 can be 
made in (^UJ, 

n = ^{a + be-">'^ + ce-^">^) . (11) 
v3 
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The emergence of cat states can be understood in terms of the information that the measurement process yields. 
Wc first consider the case of initial coherent states with well-defined phases (22| , 



(12) 



where |xP = N , i.e. the mean number of particles in each coherent state is N . For large TV, the probability density of 
atoms that would be observed if we were to take a 'snap-shot' of the overlapping clouds at a time when most atoms 
satisfy the far-field condition, is 



P{0, = ^ ( 1 + [cos(ry6l - 0ha) + cos(7/6' - ^^b) + cos(2776i - (fjba - ] , 



(13) 



where 0ba = 4'b — (f'a, 4'cb = 4>c ~ 0b, and F{ri) is the normalized probability distribution for rj. We see that, 
independent of the form of F^rj), this interference pattern is identical under exchange of (f)i,a and (pcb- This means that 
the measurement of such an interference pattern does not unambiguously establish the phase between the remaining 
condensates. Thus, the initial triple Fock state, which is an equal superposition of all relative phases evolves, after 
partial detection, into a macroscopic superposition of two sets of relative phases. This is the principle behind the 
cat generation scheme. A similar calculation for four BECs, a, b, c, and d, shows that the distribution of particles 
detected is invariant under the exchange of the phases (jjta and (f)dc- This phase ambiguity suggests that cat states 
may also be created by this technique for arrays with n = 4. It would be interesting to further investigate how cats 
emerge for n > 3. We shall now check the result for n = 3 by direct numerical calculation. 
With the first detection of an atom at 9, operator (|ll|l acts on state l(2Jl to give, 
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{\N - 1,N,N) + e-'''^\N,N - l,N) + e'^"''^\N, N,N-1)) . 



The probability distribution for the relative phases, (j)ba and (j)cb, of the remaining condensates is given by 

Pi(f>ba,(f>cb\0) = |(0ba, 0cb|V')P , 



(14) 



(15) 



where \(pba, 4'cb) is a three-mode generalization of the Pegg-Barnett relative phase state 123^ of the form, 

\4ba,M = J2 e'f>^e^'(^+*-)e''-(x+"^-+*-)b,g,r), (16) 



p,q,r 




FIG. 4: Example of the joint relative phase probability distribution for a system of three condensates (each initially in Fock 
states with 100 atoms) after 100 atoms have been detected in the far field. 
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FIG. 5: Contour plot of the number distribution of atoms in a and b corresponding to the state shown in Fig. 0] 

\p, q, r) represents a triple Fock state with p, q, and r atoms in modes a, b, and c respectively, and x can take any 
value. Substituting ((TH) and l(TB|) into l(TB|) we obtain, 

2 

P{<l>ba, (l)cb\6) = 1+3 [C0S(776' - (t>ba) + COsi-qO - (t)cb) + C0s(2r;6' - 0ba - (jicb)] ■ (17) 

As predicted, this phase distribution is symmetric under exchange of (pba and (pcb- If the distribution is peaked at 
{4'bai4>cb) = ('I'i,^2), it will also bc peaked at the converse {4>bai4>cb) — {^2,'^i)- In other words, the state is a 
superposition of these two results and so a cat state begins to emerge. This symmetry property can be shown to 
hold for subsequent detections, however the combinatorials get complicated very quickly and so it is much easier to 
numerically simulate the detection process. 

Our system for this simulation initially consists of three condensates each in a Fock state with 100 atoms. The 
condensates are partially released from the traps and detections are made on a screen in the far field. In our simulations, 
we take the detected atoms to have a range of values of rj given by the distribution, F{r]). The interference pattern 
observed on the screen has the form of 1)13(1 and so it is possible to determine the relative phases of the condensates. 
However, there is still an ambiguity as to which condensates the phases should be assigned. In Figure 4 we have 
plotted an example of the results for the relative phase probability distribution, P{(t>baT<t>cb), of the state after 100 
atoms have been detected. As expected, the state of the remaining portions of the condensates has evolved into a 
Schrodinger cat in their relative phases. The state shown here is a superposition of (jy^a ~ — 0.47r and 0cb ~ 0.27r with 
the reverse, 4>ba ~ 0.27r and 4>cb ~ — 0.47r. The positions of the peaks shown in Figure 4 differ from run to run. 

This demonstrates how it may be possible to create large Schrodinger cat states by simply measuring interference 
patterns from an array of n-tuple Fock state condensates. A particularly pleasing feature is that every trial successfully 
creates a cat state and that the interference pattern itself tells us what the two superposed phases are. This means 
that we always create a cat of a known form, but with unpredictable phases. 

During this process, there is an interplay between the detections, which tend to create relative phases between the 
condensates, and interactions, which tend to wash the phases out. It has been shown that a balance is reached between 
these two effects and so interactions between the atoms will not prevent cat states from forming. When detections 
are no longer made, the relative phases will undergo collapses and revivals _25j. The effect of the interactions could 
be reduced by making use of Feshbach resonances to tune the scattering length between atoms . 

The next challenge is to devise methods for seeing these cats directly and confirming the coherence between the 
two parts. For the system considered here, this could be achieved simply by looking for interference fringes in number 
space. We can see this by considering a state of the form depicted in Figure 4 for which (0ha, 0cfc) — (*&!, *i'2) and the 
converse. After many detections, this state can approximately be written as, 

\'4,) = 1=Y^ Cp^g (e-'P*ie*(2JV-p-?)$2 ^ g»5g-«p*2g»(2W-p-g)*i j q,2N-p- q), (18) 
P,9 
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where Cp.g is the probabihty amphtude that there are p atoms in a and q in b. There are a total of 2N particles in 
the system since N were detected in creating the state. The additional phase 6 between the terms depends on the 
particular detections made and so, although it is random, it can be known for a particular experimental run. The 
probability distribution for the number of atoms in a and b is then given by, 

P{Na,Nt) = \{Na,m,2N -Na-N,\ij)\^ 

- |C^„,jvJ'(l + cos[(27V-iVb)($i-$2)+'5]). (19) 

This predicts fringes that vary with Nb but not with Na- In Figure 5, we have plotted the number distribution for 
the cat state depicted in Figure 4 and resulting from the example simulation. This confirms that interference fringes 
of the form of (|19|l are present. In practice, this distribution would be obtained by an ensemble of measurements. We 
could create cat states similar to the one shown in Figure 4 by making measurements in the far field and, whenever a 
particular pair of relative phases is obtained, we could then directly measure the number of atoms in the remaining 
condensates. Repeating this procedure would enable us to build up the number distribution in Figure 5. No fringes are 
predicted for states that are not a superposition of different phases and so an experimental observation of interference 
would provide strong evidence that a Schrodinger cat state had been created. 
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Royal Society and Wolfson Foundation, NIST, ONR, and ARDA. 
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